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The deep X-ray, optical, and far-infrared fields that constitute GOODS are 
sensitive to obscured AGN {Nh > 10^^ cm~^) at the quasar epoch (2; ~ 2 — 3), 



^ ■ as well as to unobscured AGN as distant as z~7. Luminous X-ray emission is a 

(yQ • sign of accretion onto a supermassive black hole and thus reveals all but the most 

O 

O ■ spectral energy distributions for AGN to model the X-ray, optical and far-infrared 

flux distributions of the X-ray sources in the GOODS fields. A simple model 
based on the unified paradigm for AGN, with ~ 3 times as many obscured AGN 
2 ■ as unobscured, successfully reproduces the z-band flux distributions measured in 
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the deep HST ACS observations on the GOODS North and South fields. This 
model is also consistent with the observed spectroscopic and photometric red- 
shift distributions once selection effects are considered. The previously reported 
discrepancy between observed spectroscopic redshift distributions and the predic- 
tions of population synthesis models for the X-ray background can be explained 
by bias against the most heavily obscured AGN generated both by X-ray obser- 
vations and the identification of sources via optical spectroscopy. We predict the 
AGN number counts for Spitzer MIPS 24 //m and IRAC 3.6-8 //m observations 
in the GOODS fields, which will verify whether most AGN in the early Universe 
are obscured in the optical. Such AGN should be very bright far-infrared sources 
and include some obscured AGN missed even by X-ray observations. 

Subject headings: galaxies: active, quasars: general. X-rays: diffuse background 

1. Introduction 

Extensive studies of local Active Galactic Nuclei (AGN) have led to a unification 
paradigm wherein continuum and broad-line emission from the active nucleus are hidden 
from some lines of sight by an optically thick medium (Antonucci 1993). At such orienta- 
tions, AGN lack broad emission lines or a bright continuum and arc called Type 2 AGN (e.g., 
Seyfert 2 galaxies); usually they have strongly absorbed X-ray spectra as well. A population 
of these obscured AGN out to redshift 2-3 has been invoked to explain the X-ray "back- 
ground" (Madau, GhiseUini, & Fabian 1994; Comastri et al. 1995). Recent deep surveys 
with Chandra and XMM have resolved most or all of this background, and thus must con- 
tain high-redshift, obscured AGN (Brandt et al. 2001; Rosati et al. 2002). Obscured AGN 
are needed to explain the spectral shape of the X-ray background (Sctti & Woltjer 1989) 
since the average observed AGN spectrum (Gruber 1992) is much harder than the typical 
X-ray spectrum of an unobscurcd AGN (Mushotzky et al. 1993). Because strong absorption 
of the ultraviolet and soft X-ray emission dramatically hardens the observed spectrum of 
obscured AGN, population synthesis models involving large numbers of obscured AGN have 
been very successful at matching the X-ray background intensity and spectrum (e.g., Madau, 
GhiseUini, & Fabian 1994; Gomastri et al. 1995; GiUi et al. 1999, 2001). 

The main prediction of population synthesis models, namely that a combination of 
obscured and unobscured AGN constitute the X-ray background, has been borne out by 
deep X-ray observations (Gilli 2003; Perola et al. 2004). However, the observed redshift 
distribution of X-ray sources in deep surveys is peaked at lower redshift than these models 
require. Specifically, GiUi et al. (2001) predict a peak in the redshift distribution at z ~ 1.4, 
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and a ratio of obscured to unobscured AGN that rises from 4:1 locally to 10:1 at z ~ 
1.3. However, optical spectroscopy of X-ray sources in the Chandra Deep Fields and the 
Lockman Hole indicates a redshift peak at z ~ 0.7, and only twice as many obscured AGN 
as unobscured (Hasinger 2002; Barger et al. 2003). 

Because few Type 2 AGN are known at redshifts z ~ 2 — 3, where AGN are most 
numerous, it had been suggested that they do not exist, perhaps because the obscuring 
torus of gas and dust evaporates at high luminosity (Lawrence 1987). Now, with deep X-ray 
surveys, a few such objects have clearly been found (e.g., Norman et al. 2002; Stern et al. 
2002; Dawson et al. 2003). It is important to note that UV-excess or optical emission-line 
surveys would not have found most obscured AGN, nor would soft X-ray surveys such as the 
ROSAT All-Sky (Voges et al. 1999) or the White, Giommi & Angehni (WGA; Singh et al. 

1995) surveys. Instead, one needs to look at hard X-rays, where the absorption is smaller, 
or in the far infrared, where the absorbed energy is re-radiated. 

Discovering a previously undiscovered population of obscured AGN — a population 
suggested by the hardness of the X-ray background — was a strong motivation for the Great 
Observatories Origins Deep Survey. GOODS consists of deep imaging in the far infrared with 
the Spitzer Space Telescope (Dickinson Sz Giavahsco 2002) and in the optical with the Hub- 
ble Space Telescope (Giavalisco et al. 2004) on the footprints of the two deepest Chandra 
fields (Giacconi et al. 2001; Brandt et al. 2001; Alexander et al. 2003, hereafter A03). The 
total area is roughly 60 times larger than the original Hubble Deep Field (Williams et al. 

1996) and nearly as deep in the optical. The Great Observatories data were augmented 
with ground-based imaging and spectroscopy^. With extensive coverage over 5 decades in 
energy from 24 /xm to 8 keV (A = 1.55 A ), this survey is well suited to find a high-redshift 
population of obscured AGN if they exist. A complementary approach, given the relatively 
low surface density of AGN (compared to normal galaxies), is to target higher luminosity 
AGN over a wider area of the sky, an approach followed by, for example, the Chandra Mul- 
tiwavelength Project (ChaMP; Green et al. 2004), Calan-Yale Deep Extragalactic Research 
(CYDER; Castander et al. 2003), Serendipitous Extragalactic X-ray source identification 
program (SEXSl; Harrison et al. 2003) and the High-Energy Large-Area Survey 2 (HEL- 
LAS2XMM; Baldi et al. 2002) surveys. 

In this paper we discuss the AGN populations detected in the X-ray and optical in the 
GOODS North and South fields. Assuming a simple unification scheme, in which roughly 
three-quarters of all AGN are obscured at all redshifts, we explain the optical magnitude, 
hard X-ray flux, and redshift distributions of GOODS AGN. This model is compatible with 



Observations are summarized at http://www.stsci.edu/science/goods/ 
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previous population syntheses models for the X-ray background. We also use this model 
to predict the number counts and redshift distribution of AGN that will be detected with 
Spitzer in the GOODS fields. These predictions differ from similar calculations by Andreani 
et al. (2003) in that we include AGN evolution, which has a very strong effect, changing the 
counts by 2 orders of magnitude at the wavelengths of interest. 

In § 2 we outline the procedure used to derive the number counts and redshift distribu- 
tions at various wavelengths, and specify the AGN luminosity function and Spectral Energy 
Distributions used, which are appropriate to the unification paradigm and are based on a 
combination of observation and theory. Results arc discussed in § 3, and compared to obser- 
vations in the GOODS fields. In § 4 we present predictions for the Spitzer observations at 
24, 8 and 3.6 microns and discuss definitive tests for the obscured population. Conclusions 
are given in § 5. Throughout this paper we assume Hq — 70 km s ~^ Mpc"^, = 0.3 and 
Qa = 0.7. 



2. Calculation of Multiwavelength Number Counts 
2.1. Overview of Inputs and Procedure 

To derive the number counts at any wavelength we start with a hard X-ray luminosity 

function, an assumed cosmic evolution, and a library of spectral energy distributions. We use 
the hard X-rays as a starting point because observations at 2-10 keV in the rest frame are less 
affected by obscuration and therefore provide a less biased view of the AGN population, al- 
though they are still biased against detection of heavily absorbed sources {Nh ^ 10^^ cm~^). 
The intrinsic X-ray luminosity of each AGN is then related to its observed X-ray flux via its 
Nff value and intrinsic X-ray spectral index. 

Hard X-ray surveys are heavily dominated by AGN and thus make AGN very easy to 
identify. Although far-infrared emission can be even less biased, since the optical depth of the 
obscuring matter is low and the dust emission is quasi- isotropic (Pier & Krolik 1992), such 
surveys have a very low yield of AGN because normal galaxies are also strong far-infrared 
sources and are far more numerous. 

Hard X-ray luminosity functions based on compilations of deep Chandra, ROSAT, 
HEAO-1 and ASCA observations have been pubhshed recently by Ueda et al. (2003;U03 in 
what follows) and Steffen et al. (2003). We use the work of U03, based on 247 AGN selected 
in the hard X-ray band in deep fields hke the Lockman Hole and the Chandra Deep Field 
North. This sample covers the X-ray fiux range from 10"''^° to 3.8 x 10~^^ erg cm~^ s~^ in the 
2 — 10 keV band. The luminosity function refers to the rest-frame absorption-corrected X-ray 
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luminosity. The dependence of the luminosity function on the column density is calculated 
separately using an 'W^j function" presented in Equation 6 of U03, which is based on the 
relative number of sources at each Nh observed in their sample. 

The number of sources per unit volume per unit of logLx and per unit of log Nh is 
(U03): 

dNHCiLxdz dlogLx 
where $ is the luminosity function, which also includes evolution with redshift, / is the 
observed neutral hydrogen column density distribution and V{z) is the co-moving volume 
as a function of redshift, which depends on the adopted cosmology. We also adopt the 
luminosity-dependent density evolution model of U03, in which low-luminosity sources peak 
at lower redshift than high- luminosity AGN. This is compatible with evolution calculated in 
the optical bands by Boyle et al. (2000), which peaks at redshift 2; ~ 2 (but includes only 
high luminosity objects). We refer the reader to U03 for more details about the hard X-ray 
luminosity function. 

The Spectral Energy Distributions (SEDs) described in § 2.2 give the AGN luminosity 
at any wavelength. We use the number density in Eqn. (1) to generate a population of 
objects spanning the following ranges of Lx and Nh- Nh — — 10^^ cm~^ and Lx — 
10^^ — 10^^ ergs s"-*^. We calculate the number counts at any given wavelength by summing 
sources of the same observed flux at that wavelength, and scaling the result to the total area. 
We book keep this calculation separately for different populations, for example, unobscured 
and obscured AGN, adopting Nh = 10^^ cm~^ as the dividing point between the two classes 
(as do U03). 

The combination of the U03 luminosity function (version appropriate for our cosmology; 
sec Table 3 in U03), A'^^^ function, and the AGN SEDs described in the next section will be 
called Model A in what follows. 



2.2. AGN Spectral Energy Distributions 

The rest-frame spectrum of an AGN depends strongly on the intrinsic luminosity of the 
central engine and the amount of obscuration along the line of sight. In the simple unification 
model considered here, the obscuring matter is distributed in an axially symmetric geometry 
which, assuming random orientations of the symmetry axis, dictates the distribution of 
neutral hydrogen column density, Nh- The obscuring gas and dust both absorbs and emits 
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radiation (Nenkova et al. 2002; Elitzur et al. 2003; see also Pier & Krolik 1992, 1993; Granato 
& Danese 1994), conserving energy when integrated over all angles. 

We construct AGN SEDs from X-rays to the far infrared as a function of two parameters, 
namely the intrinsic X-ray luminosity in the 2-10 keV band and the line-of-sight column 
density of neutral hydrogen, Nh- We consider three separate wavelength regions — X-rays, 
optical/UV, and infrared — then merge the components with appropriate normalizations. 
Specifically the SEDs are constructed as follows: 

• The intrinsic relation between X-ray and UV luminosity (Vignali et al. 2003) was used 
to normalize the unobscured optical AGN spectrum, which is taken from the SDSS 
composite quasar spectrum (Vanden Berk et al. 2001). 

• Absorption was then added to both the X-ray and UV/optical parts of the spectrum. 
In X-rays photoelectric absorption was assumed, while in the UV/optical Milky- Way 
type reddening laws were used, with a standard galactic dust-to-gas ratio to convert 
Nh into optical extinction (see § 2.2.1 and 2.2.2 for details). 

• An [Mb — —20.47 mag) eUiptical host galaxy was added to the optical AGN 
spectrum. In the most obscured sources, the host galaxy dominates the rest-frame 
optical-near-IR spectrum. See § 2.2.2 for details. 

• The value of Nh for the infrared dust emission models of Nenkova et al. (2002); Ehtzur 
et al. (2003) was related to the angle between the equatorial plane of the AGN and 
the fine of sight for a simple torus geometry (§ 2.2.3). 

• The infrared spectra for different angles are normalized at 100/xm, where the emission 
from the AGN is roughly isotropic (Pier & Krohk 1992), and the infrared spectrum 
for the appropriate angle (i.e., Nh value) is added to the AGN spectrum. Details are 
presented in § 2.2.3. 

This composite spectrum is shifted to the desired redshift and the change in the SED of the 
host galaxy caused by passive stellar evolution is included (see § 2.2.2). 

Examples of the final composite SEDs for some of the X-ray luminosities and Nh val- 
ues used in this calculation are shown in Figure 1. Two known AGN, the Type 1 quasar 
PG0804+761 (^ = 0.1, A^^^ = 3.1 X 10^° cm'^; Elvis et al. 1994) and the Seyfert 2 galaxy 
NGC 7582 {z = 0.00525,A^^^ = 1.24 x 10^^ cm-^; Bassani et al. 1999) agree well with the 
appropriate model SEDs, as shown in Figure 2. Note that once Nh and redshift (and thus 
luminosity) are specified, there are no free parameters to adjust the fit to the data. 
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The effective broad-band power-law spectral slope, a^, is defined as: 



^ _ log[/.(j)//.(0] ... 

where is the rest- frame flux density and Vi and Vj are the frequencies of interest. For 
our spectral models, a source with an unabsorbed X-ray luminosity of 10^^ ergs s~^ and 
Nh — 10^° cm~^ has agx — —1.73, while a low-luminosity AGN with Lx — 10^^ ergs s~^ 
and mild obscuration, Nh ~ 10^^ cm~^. has = -1.41 (taking 2500 A and 2 keV as the 
fiducial points). Both values are well within the observed range, while the average of our 
model distribution is very similar to the average observed values (Vignali et al. 2003). 

Similarly, the values of fx/ fiR for these models, adopting the index defined by Barcons 
et al. (1995), where fx is rest-frame monochromatic absorbed fiux at 5 keV and ///j is 12- 
micron fiux, arc fx/fm = 10"^ for type 1 AGN, and fx/fm = 2.15 x 10"^ for type 2 AGN 
with unabsorbed X-ray luminosity of 5 x 10^^ ergs s~^ and = 3.2 x 10^"^ cm~^. Again, 
both indices are similar to measured values for local AGN (Barcons et al. 1995). Whether 
these SED models remain vahd at high redshift is an hypothesis that is effectively being 
tested by the comparison of predicted and observed number counts. 



2.2.1. X-ray Spectrum 
Intrinsic AGN X-ray spectra can be represented by attenuated power laws of the form: 

oc E-^e-<^^^" , (3) 

where N(E) is the number of photons with energy E; T is the power-law photon index; a{E) 
represents the cross section for photoelectric absorption of soft X-rays, given by Morrison 
& McCammon (1983) assuming solar abundance of metals; and Nh is the neutral hydrogen 
column density along the line of sight, ranging from typical high Galactic latitude values of 
Nh — 10^° cm~^ to the limit for Compton-thick absorption, Nh ~ 10^^ cm~^. A typical 
value for the intrinsic slope is F = 1.9 (e.g., Nandra & Pounds 1994; Nandra et al. 1997; 
Mainieri et al. 2002). However, it is important to note that a reflection hump can make 
the spectrum look harder, closer to F = 1.7 (Akiyama et al. 2003; Mushotzky et al. 1978; 
Nandra & Pounds 1994). Here we assume F = 1.7; we have verifled that the choice of F has 
only a minor effect in our results. 

The X-ray spectrum is normalized relative to the optical using the correlation between 
intrinsic X-ray luminosity at 2 keV and UV emission at 2500 A, Lx oc L^y (Vignah et al. 
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2003). This relation has a dispersion of ±0.06 in the exponent and is significant at the 7.9a 
level when BALQSOs are excluded (they are in any case rare) . 

2.2.2. Optical Spectrum 

Prom 1000 A to ~ 1 //m we use the Sloan Digital Sky Survey (SDSS) composite quasar 
spectrum (Vanden Berk et al. 2001), which represents an average of over 2000 SDSS quasars 
with a median redshift z — 1.253, covering rest-frame wavelengths from 800 A to 8555 A at 
a resolution of ~ 1 A . The intrinsic luminosity of quasars in this sample spans the range from 
Mr' = — 18 mag to M^/ = —26.5 mag. This spectrum well represents unobscurcd (type 1) 
AGN, in which the optical light from the central engine is not absorbed by the dust or gas 
torus. Given that the SDSS quasar composite includes AGN brighter than the average AGN 
observed in X-rays, the use of the SDSS composite spectrum may not be realistic; however, 
it is completely appropriate to our very simple model which assumes the obscuring torus 
is independent of luminosity or redshift. Furthermore, there are good indications that the 
optical spectrum of faint quasars is very similar to the spectrum of the average SDSS quasar 
(Steidel et al. 2002). 

To calculate the effects of absorption on the optical and UV spectrum we use the ex- 
tinction coefficients of Cardelli et al. (1989), obtained from absorption measurements on 
Milky Way stars. A standard Milky Way dust-to-gas ratio was assumed to convert neutral 
hydrogen column density into optical extinction: Nh — 1.96 x lO^^Ay (Bohlin et al. 1978). 
Some studies of the dust-to-gas ratio in AGN point out that this ratio is on average signif- 
icantly smaller than the value obtained for the Milky Way (Maiolino et al. 200 lb, a), but 
this picture has been questioned (Weingartncr & Murray 2002). Here we assume the more 
conservative position and use the standard dust-to-gas ratio value. If in fact this ratio is 
two orders of magnitude lower for AGN (Maiolino et al. 2001b), marginally obscured AGN 
(with Nff ~ 10^^ cm~^) will be significantly brighter in the optical. For example, a source 
with Lx — 10^^ ergs s"^ and Nh — 10^^ cm~^ at 2; = 1 will be ~1 magnitude brighter in the 
2;-band if the dust-to-gas ratio is two orders of magnitude smaller than in the Milky Way. 
For more luminous sources this effect is larger, reaching 5 magnitudes in the 2;-band for a 
source with Lx = 10"^^ ergs and Nh = 10^^ cm^^ at z = 1. However, this is about the 
maximum discrepancy for the sample studied in this paper, since sources with column densi- 
ties larger than ~ 5 x 10^^ cm^^ or lower than 10^^ cm^^ arc not significantly affected in the 
optical bands since they arc already too obscured (and thus dominated by the host galaxy) 
or the obscuration is too low to make a significant change. Also, the infrared bands are not 
greatly affected by a change in the dust-to-gas ratio assumed. Therefore, our conclusions 
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should not be greatly affected by a different choice of dust-to-gas ratio. 

Given the small GOODS volume, the contribution of very high redshift sources {z > 5) 
to the total sample is negligible. For these sources, z-band observations (Agg^ ~ 8800 A ) 
effectively correspond to A > 1460A in the rest frame, and therefore absorption from the 
intergalactic medium (Madau 1995) can safely be ignored in the whole sample. 

To the AGN spectrum we add optical hght from a host galaxy, for which we assume an 
elliptical with spectrum given by Fioc & Rocca-Volmerange (1997). In order to account for 
the evolution of the host galaxy, caused mainly by star formation, we adopt the evolutionary 
correction of Poggianti (1997), assuming an elliptical galaxy with an e-folding time of the 
star formation rate of 1 Gyr. This model predicts an increase in the 2;-band himinosity of 
~ 1.2 magnitudes at redshift ~ 1. Preliminary results using GOODS optical data (Simmons 
et al 2004, in prep) show that a large fraction of the AGN host galaxies in the survey are 
luminous elliptical, justifying our assumption. The host galaxy dominates the optical light 
for obscured sources with Nh > 10^^ cm~^. As an example, for a source with intrinsic hard 
X-ray luminosity Lx — 10^^ ergs s~^ and Nh — 10^^ cm~^ the AGN contribution to the total 
optical hght is only ~ 4% a,t z — 0. 



2.2.3. Infrared Spectrum 

Optical, UV, and soft X-ray light absorbed by dust within the AGN is re-radiated as 
thermal emission in the far infrared. Instead of using observed infrared spectra we use the 
dust- re-emission models of Nenkova et al. (2002); Ehtzur et al. (2003) in order to construct 
a grid of infrared spectra as a function of the X-ray luminosity of the AGN and observed 
Nh value (related to the hne-of-sight angle in these physically-motivated models). These 
models provide a better fit for recent infrared observations, as well as other advantages over 
older torus models (e.g. Pier & Krolik 1992, 1993; Granato & Danese 1994); for a detailed 
discussion of this point see Elitzur et al. (2003). 

The Nenkova et al. (2002) models postulate a random distribution of clumps inside a 
dusty torus. Each clump is optically thick, and the radial distribution is confined between an 
inner radius Ri and an outer radius Ro (see Figure 2 of Elitzur et al. 2003). The mean free 
path between clumps is given by a power law of the form r"^ where r is the radial distance 
from the central engine, while the angular dependence of the number of clumps is a Gaussian 
distribution with dispersion cr, so that the number of clouds as a function of viewing angle 
(5 is given by 

iVr(/3) = Arr(0)exp('-^') , (4) 



-10- 



where (5 is measured with respect to the equatorial plane. In this work we assume Nt{0) — 10, 
(7 = 29° and an optical depth of each clump of ry = 100, which gives an optical depth at the 
equator of 1000 (or = 1-2 x 10^^ cm^^ assuming the standard dust to gas ratio) and 0.06 
at the poles. This permits AGN ranging from those that are Compton thick down to those 
that are completely unobscured (although to compare to observations, we assume a minimum 
column density of Nh ~ 10^° cm~^ for unobscured sources, corresponding to the observed 
column density through our own Galaxy in the direction of the GOODS-S or GOODS-N 
fields). Also, we consider three different combinations of values for the parameters Ri/Rp 
and q; namely Ri/Ro — 30 and q — I, Ri/Ro — 30 and q — 2 and Ri/Ro — 100 and q — 1 
(Elitzur et al. 2003). All three combinations of values for the parameters generate IR spectra 
that are consistent with the observations of some local AGN, as can be seen in Figure 2 and 
in Nenkova et al. (2002). This dust emission model gives an IR spectrum for a given intrinsic 



We calculate the Nh distribution expected from a simple obscuring torus model with 
fixed geometry and dust distribution and we use this model to convert Nh into viewing an- 
gle. From this model we obtain a new Nh function, based on the AGN unification paradigm, 
called model B in what follows (in contrast to model A, which uses the observed Nh distri- 
bution). We assume the torus lies at a distance Rm from the central engine and has height 
Rt, and that its density distribution is given by 



(no radial dependence). A schematic diagram of this geometry is shown in Figure 4 and the 
optical depth as a function of viewing angle is 



where Te is the optical depth of the torus at the equatorial plane and 7 parameterizes the 
exponential decay of density with viewing angle. Using a standard dust-to-gas ratio, we can 
relate optical depth to column density via Nh — 1.2 x lO^^r cm~^ (Bohlin et al. 1978). Then 
for an equatorial column density Nh — 10^^ cm~^, a ratio of radius Rm/ Rt — 1-01, 7 = 4, 
and a random distribution of viewing angles, we obtain the model Nh distribution shown in 
Figure 3. 

This set of parameters was selected in order to obtain a ratio of obscured to unobscured 
of ~ 3, consistent with the locally observed value (Risahti, Maiolino, & Salvati 1999) and 
roughly consistent with population synthesis models that explain the spectrum of the X-ray 
background. A larger torus at larger radius could also give a 3:1 ratio and would have a cooler 



Nh. 



p{9) oc exp(— 7I cos^l) , 



(5) 



t{9) 
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far-infrared spectrum, but until we get the Spitzer infrared data, we cannot constrain the 
far-infrared emission. This distribution is very similar to that observed in the GOODS fields 
(solid line in Figure 3), except for > 10^'^ cm^^ sources. For such high column densities, 
the discrepancy, of the order of 10-15% in the fractional distribution, is probably caused by 
incompleteness in the Chandra samples, since the amount of obscuration is large enough to 
hide even the hard X-ray emission. In particular, for Nh — 3 x 10^^ cm~^ the incompleteness 
of the Chandra observations in the Chandra Deep Fields is ~ 25%, increasing to ~ 70% for 
Nff > 10^^ cm~^, using the U03 luminosity function and calculating the number of sources 
below the X-ray flux hmit as a function of Nh- Overall, for our model, roughly 50% of the 
AGN are not detected in the Chandra Deep Fields (see § 4.3 for details). 

Model B thus comprises the U03 luminosity function (including its redshift evolution), 
the AGN SEDs described in section §2.2, and the Nh function given by equation 6. Implicit 
in this model is that the ratio of obscured to unobscured AGN is constant with redshift, and 
that the geometry of the torus does not change with redshift or luminosity. 

For a given luminosity, spectra for different orientation angles are normalized at 100 /im, 
where the optical depth is low and thus the re-processed emission from the torus is roughly 
isotropic. The infrared spectrum is then joined smoothly to the optical spectrum at 1 /im 
so that the resulting spectrum is continuous. Infrared emission from star formation is not 
considered in this model. 

3. Observed and Predicted AGN Number Counts 

3.1. The GOODS Data 

We compare the number counts calculated above to the optical and X-ray flux distri- 
butions for AGN in the two GOODS flelds, based on Chandra and HST ACS multi-band 
data. The two GOODS fields, each 10' x 16', were imaged with ACS in the B,V,i and z 
bands, for a total of 13 orbits per pointing, reaching AB magnitude 27.4 (5(t) in the z band^. 
The GOODS fields have also been imaged extensively from the ground in UBVRIzJHK, and 
spectroscopy is ongoing. Details of these observations can be found in Giavalisco et al. (2004) 
and at http://www.stsci.edu/science/goods. 

GOODS-N and GOODS-S fields have pubhshed deep X-ray observations: the 2 Ms 



^Observations taken with the NASA/ESA Hubble Space Telescope (HST), which is operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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Chandra Deep Field-North (A03) and the 1 Ms Chandra Deep Field-South (Giacconi et al. 
2002, hereafter G02). These two ultra-deep X-ray surveys provide the deepest views of the 
Universe in the 0.5-8.0 keV band. The CDF-N is ~ 2 times more sensitive at the aim 
point than the CDF-S, with 0.5-2.0 keV and 2-8 keV flux limits (S/N= 3) of ^ 2.5 x 
10"^'' erg cm^^ s^^ and ^ 1.4 x 10^^^ erg cm~^ s~^, respectively (AOS). The CDF-N remains 
> 1.8 times and > 1.5 times deeper than the CDF-S over ?a 50% (~ 90 arcmin^) and ~ 75% 
(Ri 135 arcmin^) of the area of GOODS fields, respectively. Point-source Chandra catalogs 
have been produced by AOS for the CDF-N and CDF-S and by G02 for the CDF-S. In this 
study we use the Chandra catalogs of AOS which were generated using the same methods for 
both fields. These catalogs include S26 sources detected in the CDF-S and 50S in the CDF- 
N. 22S of the X-ray sources in the CDF-S are located in the GOODS-S region, while S24 of 
the sources in the CDF-N were found in the GOODS-N region. If only the sources detected 
in the hard band are included, the sample is reduced to 141 sources in the GOODS-S field 
and 210 sources in the GOODS-N region. 

The GOODS fields will be observed with the Spitzer IRAC instrument in all four bands 
(S.6 to 8 /im), with expected sensitivity of 0.6/iJy in the S.6/im band. Both fields will also 
be observed with the Spitzer MIPS instrument at 24 microns, to a fiux limit that depends 
somewhat on the as- yet unknown source density (and hence confusion limit), but which we 
take to be roughly 22 /xJy {5a). 

There are 168 published spectroscopic redshifts for X-ray sources in the CDF-S, all of 
them measured using the FORSl and F0RS2 cameras at the VLT telescopes (Szokoly et al. 
2004). Photometric redshifts have been calculated for all the sources detected in the optical 
bands in the GOODS South field (see Mobasher et al. 2004 for details), which account 
for ~ 90% of the observed X-ray sources. Properties of some of the remaining sources 
not detected in the ACS observations are described in detail by Koekemoer et al. (2004). 
Spectroscopic redshifts were used when available. 

Spectroscopic redshifts for 284 X-ray sources in the CDF-N were obtained by Barger 
et al. (200S) using the Low-Resolution Imaging Spectrograph and the Deep Extragalactic 
Imaging Multi-Object Spectrograph on the Keck 10-m telescopes and the HYDRA spectro- 
graph on the WIYN S.5-m telescope. Photometric redshifts were used when spectroscopic 
redshifts were not available (see Barger et al. 200S). Spectroscopic redshifts were obtained 
for 54% of the sources, while if photometric redshifts are added the completeness level rises to 
~ 75%. Sources without spectroscopic or photometric redshifts are very faint in the optical, 
with all but two having z > 24.0 mag. Given the very faint magnitudes of the sources with- 
out redshifts in the GOODS-N region, they are likely at redshift z > 1.0 and lack the strong 
broad emission lines normally used to identify unobscured AGN at high redshift. That is. 
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most are probably obscured AGN at z > 1 (e.g., Alexander et al. 2001; Koekemoer et al. 
2002). 

3.2. The AGN Sample 

Almost all the X-ray sources detected in the GOODS fields are AGN, with the exception 
of a few nearby starburst galaxies (Alexander et al. 2003). Sources in the X-ray catalogs 
were matched to the ACS images of the GOODS North and South fields by Bauer et al. 
(2004, in prep) using the likelihood method described in Bauer et al. (2000). Using only 
X-ray sources that were detected in the hard band and were unambiguously identified in the 
optical images (i.e., only one optical counterpart within ~ 1 arcsecond of the X-ray centroid), 
the final sample includes 128 sources in the GOODS South region and 178 in the GOODS 
North field. It is important to note that reducing the analysis to the unambiguously matched 
sources eliminates ~ 10% of the hard X-ray sources in the GOODS regions. Roughly half of 
these sources in the GOODS-S field have counterparts in deep X-band imaging and thus may 
be obscured AGN at high redshift (Koekemoer et al. 2004). Optically faint X-ray sources 
detected in near-infrared bands in other fields were also discussed by Mainieri et al. (2002) 
and Mignoh et al. (2004). 

3.3. Deriving the Nh Distribution from the Chandra Data 

The hardness ratio, which we calculate for each source, is defined as {H — S)/{H + S), 
where S is the number of counts detected in the 0.5 — 2 keV band and H is the number of 
counts in the 2 — 8 keV band. The neutral hydrogen column density was calculated from 
the hardness ratio assuming an intrinsic power-law X-ray spectrum with photon indices 
r = 1.7 or r = 1.9. We generated a conversion table using XSPEC (Arnaud 1996) to 
calculate hardness ratios for a range of Nh and redshift values; this program incorporates 
the Chandra ACIS instrumental response matrix. We added Galactic absorption column 
densities (Stark et al. 1992) at z = of A^if = 8 x lO^'^ cm-^ for the GOODS-S field and 
Nh — 1.6 X 10^° cm~^ in the GOODS-N field to each Nh-z pair stored in the table, since 
all X-rays pass through the interstellar medium of our galaxy. Using the table, the hardness 
ratio of each source with a spectroscopic redshift (photometric redshifts are too uncertain 
for this purpose) can be converted to a value of Nh- The distribution of derived Nh values 
for 82 sources with spectroscopic redshifts in the GOODS-S field and 103 sources in the 
GOODS-N field is shown in Figure 3. Note that for individual sources, our method of 
deriving Nh may not be very accurate, as some will have soft excess emission or intrinsically 
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soft spectra (indeed, these may be the reasons why there is an excess of apparently low 
column density objects) or will have intrinsically steeper spectra (and thus have spuriously 
high Nh) and therefore our resulting distribution can be affected by systematic errors. In 
order to minimize these effects, the value of Nh determined for each source is only used 
to calculate the distribution of the sample and not for other purposes (e.g. to correct the 
observed X-ray flux). 

The predicted Nh distributions for the sources in the GOODS N+S fields in models 
A and B arc compared in Figure 3. Model B provides a better fit to the observed data; 
a K-S test on both models compared to the observations revealed that the null hypothesis 
(that the model and observed distributions are drawn from the same parent distribution) is 
acceptable at the 75.3% confidence level for model B and at the 58.5% confidence level for 
model A. This is not surprising since model B uses an Nh function with parameters chosen 
to be consistent with the observations; however, it is important to note that this model 
is motivated by the unification paradigm, which clearly can account for the observed Nh 
distribution. 



3.4. X-ray Number Counts 

The total area covered as a function of X-ray fiux hmit for each GOODS field was 
calculated based on the results presented in A03 and is shown in Figure 5. These curves 
are needed to normalize the hard X-ray flux distribiitions for sources in the GOODS fields, 
shown in Figure 6. We also plot the X-ray distributions calculated from the hard X-ray 
luminosity function U03 for both models A (U03 observed Nh function) and B (our Nh 
function based on a simple unified model). The agreement is very good for both models; 
a K-S test shows that the null hypothesis is acceptable at > 90% confidence when either 
model was tested against the observed distribution. Although the U03 sample includes the 
CDF-N, this is not a circular argument since in that work only the brightest CDF-N sources 
(with /x > 3 X 10~^^ ergs cm~^ s~^ in the hard band) were considered, and therefore our 
calculation tests the extrapolation to much fainter fiuxes. In the hard X-ray fiux distribution 
both unobscured (type 1) and obscured (type 2) AGN are well represented (dotted and 
dashed lines in Figure 6). 
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3.5. Optical Number Counts 

In deep optical imaging with the HST ACS camera the vast majority of the X-ray 
sources in the GOODS fields are detected. Figure 7 shows the distribution of observed z- 
band magnitudes for the X-ray sources in the two fields combined (solid line), along with 
the model predictions (dashed line). The agreement between the observed and predicted 
distributions is very good; a K-S test comparing the observed distribution to the prediction 
of model B gives a confidence level for the null hypothesis of 71%, while for model A the 
confidence level is 53%. The observed distribution is very broad, with the brightest objects 
at ^ ~ 17 mag and the faintest below ^ ~ 28 mag. Unobscured AGN (dashed line) fail to 
account for the faintest optical counterparts. Given their large X-ray to optical flux ratios, 
the fainter optical sources [z > 23.5 mag) may be obscured AGN at redshifts 2; ~ 1 — 3 
(Alexander et al. 2001; Koekemoer et al. 2002). Indeed in our models most obscured AGN 
(dotted fine) have faint optical magnitudes, while the unobscured AGN are responsible for 
the peak at brighter magnitudes. For z > 23.5 mag, the approximate limit for ground-based 
spectroscopy, obscured AGN are the dominant population. 

3.6. Redshift Distribution 

Figure 8 shows the redshift distributions for the sources in the GOODS-North and South 
fields (thick solid fines) compared to the expected redshift distributions from our model B 
(dashed fines). Using photometric redshifts (Mobasher et al. 2004), our GOODS-S sample of 
hard X-ray sources is 100% complete. However, in the GOODS-N field, combining spectro- 
scopic and photometric redshifts (Barger et al. 2003), the sample is only ~ 75% complete^. 
Spectroscopic redshifts are shown by the hatched regions. At redshifts above z ~ 1, there is 
a clear discrepancy between the spectroscopic redshifts and either the photometric redshift 
or the model predictions, in the sense that there are fewer AGN with high spectroscopic 
redshifts. This is explained at least in part by the effective brightness limit for spectroscopy, 
i? < 24 mag, for even the largest ground-based telescopes. Imposing this optical limit on 
the GOODS AGN model (long dashed fine in Figure 8) does match the observed spectro- 
scopic redshift distribution very well. Obscured AGN at 2; > 1, in particular, are fainter 
than this limit and thus are not included in the spectroscopic samples. This likely accounts 
for the discrepancy between the observed redshift distribution of X-ray sources (Hasinger 
2002; Szokoly et al. 2004) and the distribution predicted from population synthesis models 
for the X-ray background (Gilfi et al. 1999, 2001), as well as the low number of Type 2 AGN 



^We expect the photometric redshifts to be 100% complete once the GOODS ACS data are fully analyzed. 
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identifications. 

The agreement between the model and the observations is good, with a K-S test giving 
63% confidence for the null hypothesis in the GOODS-S field and 21% confidence in the 
GOODS-N field. This shghtly lower level arises because the model predicts still more high 
redshift AGN than are observed in the GOODS fields. This is because some high-redshift 

obscured AGN arc too faint even for the HST images (see Koekemoer et al. 2004) and because 
the most obscured AGN are not detected in the Chandra deep fields. Also, an excess of 
observed sources sX z < 1 can be seen. This is explained by the presence of clusters and 
large scale structure in both fields (e.g. Gilli et al. 2003). 

Figure 8 shows the difference between the North and South fields, in the sense that 
there is a larger discrepancy between the photometric redshifts and the model in the North. 
This must be due to incompleteness since again the missing 25% of the AGN (those without 
photometric or spectroscopic redshifts) are preferentially fainter. 

The model redshift distribution is therefore in agreement with the data once these 
selection effects are considered. It is also compatible with the kind of distribution predicted 
from population synthesis (Gilli et al. 2001), in the sense that it peaks at a higher redshift 
than previously reported distributions, with the same caveat about selection effects. This 
better agreement between models and observations is explained also in part by the use of the 
U03 luminosity function, that includes a redshift distribution that peaks at lower redshifts 
than existing optical quasar luminosity functions (e.g. Boyle et al. 2000). 

4. Discussion 

Hard X-ray surveys find obscured sources that are largely missed in deep optical/UV 
surveys. Deriving a hard X-ray luminosity function or a redshift distribution imposes an 
effective optical cut at i? < 24 mag, since optical spectroscopy is required to obtain redshifts. 
Thus published redshift distributions can be missing optically faint sources. These sources 
can be either low-luminosity unobscured AGN or obscured AGN, depending on their central 
engine luminosity, and in general, the spectroscopic incompleteness increases with redshift. 
This selection effect explains the discrepancy between the photometric and spectroscopic 
redshift distributions (Fig. 8). 

Evidence for the existence of a significant number of obsciued AGN at ^ > 1 was 
previously given by Fiore et al. (2003), who found a correlation between X-ray luminosity 
and X-ray-to-optical flux ratio in obscured AGN. Such a correlation can be explained by our 
models since for obscured AGN the optical light is dominated by the host galaxy, which is 
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independent of the AGN luminosity. Therefore, for an obscured AGN the optical emission is 
roughly constant, while the X-ray emission scales directly with the AGN luminosity. Using 
this correlation, Fiore et al. (2003) estimated source redshifts using just the X-ray-to-optical 
flux ratio, thus including sources too faint for optical spectroscopy. They concluded there is 
a significant number of obscured AGN aX z > 1, and that these sources will be missed by 
surveys that rely on optical spectroscopy. We arrive at the same conclusion from a different 
direction, based on the comparison of our model with multiwavelength observations of X-ray 
sources. 



4.1. Optically Faint X-ray Sources 

The GOODS HST data reveal an appreciable number of optically faint, X-ray-sclccted 
AGN. While unobscured sources are bright in the optical bands, the vast majority of them 
brighter than the spectroscopic limit, obscured sources will be optically fainter and therefore 
missed preferentially by surveys that depend on spectroscopic identifications. Most obscured 
AGN are bright enough in hard X-rays to be detected in the Chandra observations. There- 
fore, given the large X-ray-to-optical fiux ratios of the optically faint sample, their hardness 
ratio and their red colors, they are very likely to be obscured AGN at z > 1 (Alexander et al. 
2001; Koekemoer et al. 2002). Sources with high X-ray-to-optical fiux ratio in other fields 
were discussed previously in the literature (e.g., Fiore et al. 2003; Mignoli et al. 2004; Gandhi, 
Crawford, Fabian, & Johnstone 2004), and in most cases these sources can be identified as 
obscured AGN at high redshift. 



4.2. Nh Distribution 

Both models for the Nh distribution fit the observed fiux distributions well (Figs 6,7). 
Model A assumed the empirical Nh function described in U03, which is based on observations 
from several X-ray surveys and did not go as deep as the GOODS sample, while model B uses 
an Nh function based on a simple unified model, with the idea of a dust torus covering the 
central engine. It is important to note that model B predicts a larger number of sources with 
Nh ~ 10^^ cm~^ than is observed, commensurate with the large number of very obscured 
sources needed to explain the X-ray background spectrum (Gilli et al. 2001, U03). The 
observed ratio of obscured to unobscured sources in the GOODS fields is ~2.5 when the 
division point is set at Nh = 10^^ cm^^ as assumed by Gilli et al. (2001), less than the 
intrinsic value of 4 required by the X-ray background models. That is, the observed ratio is 
affected by incompleteness in the X-ray samples at high column densities. Based on the U03 
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luminosity function and our torus model to calculate the Nh distribution we calculate that 

the completeness level of the Chandra observations in the GOODS fields drops to ~ 75% 
for Nh = 10^^ cm~^ and ~ 30% for Nu = 10^^ cm~^. This is caused by absorption of the 
X-ray emission which makes harder to detect sources with higher column densities in a fiux- 
limited survey. Therefore, if there are heavily obscured AGN in the field, the intrinsic ratio 
of obscured to unobscured sources is larger than the observed ratio, and can even be ~ 4, as 
suggested by population synthesis models for the X-ray background. The relation between 
the observed X-ray sources in the GOODS fields and the models for the X-ray background 
will be analyzed in more detail in a later paper (Treister et al. 2004, in prep.). 

We can ask whether the X-ray absorption correlates with the optical dimming; that 
is, is obscuration making bright hard X-ray sources optically faint? The relation between 
optical magnitude and amount of obscuration is shown in Figure 9. Sources that present 
large amounts of obscuration in the X-ray spectrum are in general the fainter sources in the 
optical, as expected. There is no strong correlation present in Figure 9 — both obscured 
and unobscured AGN appear at bright magnitudes — but the highest column densities do 
correspond to the faintest magnitudes. In particular, it is important to note that there are no 
unobscured sources at magnitudes i > 24.5, which is consistent with our model predictions 
of a very low number of unobscured AGN at faint magnitudes. There is no strong correlation 
because the optical emission from obscured AGN is dominated by the host galaxy; once the 
optical emission of the AGN is absorbed, the i-band magnitude becomes insensitive to the 
amount of obscuration or to the intrinsic luminosity of the AGN. 

The full X-ray (0.5-8 keV) to optical {z-ha.nd) fiux ratio"^ for GOODS AGN is shown in 
Figure 10, together with the ratio obtained from our model. The agreement between these 
two distributions is very good, with a K-S confidence level of ~ 92% for the null hypothesis, 
except for a small discrepancy at logF^/Fopt > 0.5, which is consistent with the previously 
calculated incompleteness of the X-ray samples for > 10^^ cm~^. A second discrepancy 
appears at the low Fx/^opt end, which can be explained by the presence of a few starburst 
galaxy and other non-AGN X-ray sources in the Chandra Deep Fields observations (e.g., 
Hornschemeier et al. 2001; Alexander et al. 2002; Bauer et al. 2002). These sources, around 
10% of the total sample, are not accounted in our model and therefore increase the number 
of observed X-ray sources with logF^/Fopt < — 2. This number of non-AGN X-ray sources 
is in agreement with the values reported by A03. 



^Defined as log Fx /i^opt = log Fx + 0.4 x -|- 4.934, where Fx is the X-ray flux in units of ergs cm ^ s ^ 
and ruz is the optical magnitude in the z band in the AB system. 
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4.3. Predictions for the Infrcired Number Counts 

The Spitzer Space Telescope was launched in August 2003 and the GOODS fields 
are scheduled to be observed in 2004. AGN are luminous infrared sources (Sanders et al. 
1989), so the ratio of obscured to unobscured AGN will be strongly constrained by the 
observed Spitzer number counts of X-ray sources. These two classes can be separated using 
a combination of IRAC 4.5 and 8 microns bands and MIPS 24 microns data (Andreani et al. 
2003). Furthermore, the infrared spectrum is sensitive to parameters of the dust emission 
model, so will help refine the present simple spectral models. 

The AGN number counts at 24 //m calculated according to § 2 assuming the appropriate 
hard X-ray fiux limit for each Chandra deep field separately, and normalized to the total 
GOODS area (0.08 deg^), are shown in Figure 11; the number counts at 8 /xm and 3.6 iim 

are presented in Figures 12 and 13. Contributions from obscured and unobscured AGN 
were calculated from hard X-ray luminosity functions using model B. For the infrared part 
of the SED, we assumed three difi^erent sets of parameters, as is described in §2.2; each 
shaded region for the total and obscured number of sources corresponds to the extremes 
that these values can have. For unobscured sources, the shape of the infrared SED is almost 
independent of the assumed parameters, so only one line is plotted. 

One important conclusion from this calculation is that all the AGN detected in the 
Chandra deep fields should be detected in the Spitzer observations. The contrary is not 
the case since some obscured AGN are missed by X-ray observations. This difference is 
shown by the top two curves in Figure 11. Indeed, obscured AGN with Nh > 10^^ cm^^ 
should all be detected with Spitzer, in principle allowing for a complete study of the AGN 
population in this field and providing a test for both the unified model of AGN and the 
population synthesis models used to explain the X-ray background. However, to identify 
obscured AGN not detected in X-rays and to separate them from luminous starburst galaxies 
will not be easy. In fact, Andreani et al. (2003) reported that is not possible to separate 
very obscured AGN activity from a burst of star formation just on the basis of infrared 
photometry, even if information in the MIPS 70-micron band is included. This problem 
is even more serious since typically there is overlap between the two populations, and a 
starbust galaxy can harbor a very obscured AGN. In order to solve this problem, we can 
take advantage of GOODS multiwavelength coverage, ranging from 0.4 to 24 microns, which 
allows us to calculate accurate photometric redshifts and bolometric luminosities for even 
very faint sources. Then, to first order, we can use bolometric luminosity in order to separate 
AGN activity from starburst galaxies, since all AGN should have L^qI > 10^^ ergs s~^, while 
starbursts typically have much lower luminosities. 

Similar calculations of the AGN number counts in the infrared were performed recently 
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by Andreani et al. (2003). For Qa = 0.7 and Qrn = 0.3, they found a total of ~ 3 x 10^ 
sources per square degree with a 24 /im flux higher than 550 /iJy, ~ 10% of them Seyfert 2s. 
We predict a much lower number of sources, ~ 10^ sources per square degree at the same 
flux limit in the 24 fim band, with ~ 70% of them being Seyfert 2s. This large discrepancy 
of 2 orders of magnitude can be explained by the difference in the assumed luminosity 
function and evolution. We used a hard X-ray luminosity function and luminosity- dependent 
density evolution , which peaks at redshifts 2; ~ 2 for QSO-like objects and 2; < 1 for low- 
luminosity AGN, whereas Andreani et al. (2003) used a local luminosity function based on 
12-//m observations and optical QSO evolution with a peak at 2; ~ 3 and exponential decay 
until z — 10. Locally both luminosity functions are very similar so the major part of the 
discrepancy is explained by the difference in the assumed evolution. This shows clearly 
that Spitzer observations will provide a significant constraint on the AGN number density 
evolution. 



5. Conclusions 

In this work we modeled the AGN population of the GOODS fields in order to explain 
the observed numbers and brightnesses at optical and X-ray energies. We also predicted 
the number counts at infrared wavelengths that will be measured with the Spitzer Space 
Telescope. Basic ingredients of our models are few. First, we used the hard X-ray luminosity 
function and luminosity-dependent evolution determined by U03. Second, we developed a 
library of multi-wavelength AGN SEDs parameterized as a function of only two parameters, 
the intrinsic rest-frame hard X-ray luminosity and the amount of gas and dust in the line 
of sight, given as a neutral hydrogen column density. These composite SEDs assume an 
intrinsic power law plus photoelectric absorption in X-rays; a QSO composite spectrum in 
the optical, independent of both redshift and luminosity; Milky Way-type extinction; an 
elliptical host galaxy evolving with redshift; and infrared dust re-emission models from 
Nenkova et al. (2002) and Elitzur et al. (2003). Finally, we needed to assume the relative 
number of sources at each value of Nh- In model A, we used the Nh function given in U03, 
which is based on observations of X-ray sources in several surveys. For model B, we used 
a simple geometric model for a dust torus to generate the expected distribution of Nh for 
random orientation angles. The torus parameters were chosen so that the ratio of obscured 
{Nff > 10^^ cm~^) to unobscured AGN is 3:1, and we assumed it was independent of redshift 
and luminosity. With these ingredients, we calculated the expected distribution of sources 
at different wavebands and compared them to the GOODS observations. 

We found general agreement between the models and observations, especially for the very 
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simple unified model. This means there may well be a significant population of obscured AGN 
at high redshift, consistent with a roughly constant ratio of obscured to unobscured AGN out 
to high redshift, z ~ 2 — 3. These objects will be bright enough to detect with the planned 
Spitzer observations in the GOODS field. The agreement between model and observations is 
remarkable given the very simple assumptions used in the calculation. Because the GOODS 
observations are so deep, there are large extrapolations from existing luminosity functions. 
The excellent agreement supports the unified model of AGN. 

A large population of obscured AGN in the early Universe explains not only the GOODS 
data but the X-ray background spectral shape and intensity. However, for this picture to 
be correct, existing X-ray samples, even from the extremely deep Chandra fields, must be 
incomplete. We find that the luminosity-function-weighted fraction of obscured AGN that 
is missed at current detection limits is approximately equal to the excess predicted by our 
unified model at faint optical magnitudes and high absorbing column densities. Spitzer 
observations, together with follow-up studies of host galaxy morphologies and environments, 
offer the opportunity to test and refine this picture. 
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Fig. 1. — UV to far IR spectral energy distributions for AGN with three intrinsic hard X-ray (2-8 keV) 
himinosities, Lx, and neutral hydrogen column densities, Nh, ranging from 10^*^ cm~^ {light lines) to 
10^^ cm~^ {heavy lines). Details about these model SEDs can be found in § 2.1. 
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Fig. 2. — Spectral Energy Distributions for the Type 1 quasar PG0804+761 (upper crosses) 
and the Seyfert 2 galaxy NGC 7582 (lower crosses). Over- plotted are our model SEDs 
(solid and dashed lines) for the appropriate Lx and Nh', an X-ray luminosity of Lx = 
8.5 X 10^5 ergs s"^ and Nh = 3 x lO^^ cm-^ for PG0804+761 (Elvis et al. 1994) and 
Lx = 9.4 X lO'^i ergs s'^ and iV^^ = 1.2 x 10^^ cm^^ for NGC 7582 (Bassani et al. 1999). 
With no free parameters to adjust, the agreement between the model SED and the observed 
values is remarkably good. 
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Fig. 3. — Distribution of intrinsic obscuring column density for X-ray sources in the GOODS fields 
with spectroscopic redshifts. The Nh values were derived assuming an intrinsic power-law spectrum 
with photon index F = 1.9 {solid line) or F = 1.7 {dashed line) and photoelectric absorption cross 
sections given by Morrison Sz McCammon (1983), shifted to the rest frame using redshifts from 
(Szokoly et al. 2004). Dashed-dotted line: Nh distribution of U03 (model A). Light solid line: 
Nh distribution calculated assuming a simple unified model with a dust torus geometry for the 
obscuring material and the GOODS flux limit (model B). As can be seen in this plot, model B 
(which is physically motivated) provides a good fit to the observed sources in the GOODS fields, 
apart from the expected incompleteness at high column densities. 
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Fig. 4. — Schematic diagram of the AGN torus model used to calculate the intrinsic Nh 
distribution expected for a random orientation. In this plot Rm is the distance from the 
black hole to the center of the torus and Rt is the radius of obscuring material. The torus 
parameters are chosen so that equatorial lines of sight correspond to column densities of 
Nh — 10^^ cm~^ and polar hues of sight have obscuration comparable to or less than typical 
Galactic values, Nh ~ 10^° cm~^. 
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Fig. 5. — Area covered as a function of depth reached in the GOODS North (2 Ms) and 
South (1 Ms) fields with the Chandra X-ray observatory, based on the calculation presented 
in A03 over only the GOODS regions. This covered area as a function of depth was used in 
the subsequent calculation of expected number counts. 
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Fig. 6. — Hard X-ray flux (2-8.0 keV) distribution for the entire sample of hard X-ray 
detected sources in the GOODS-North and -South fields {heavy solid line, sources with or 
without spectroscopic redshifts), compared to the number counts calculated from the hard 
X-ray luminosity function and models A {dot-dashed line) and B {light solid line), with the 
individual contributions of unobscured (Type 1) {dashed line) and obscured (Type 2) {dotted 
line) AGN shown separately. Note that Type 1 and Type 2 AGN follow similar distributions, 
as expected from the unified model, since hard X-rays are not strongly affected by absorption. 
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Fig. 7. — Distribution of observed z-hand magnitudes for the entire sample of GOODS-North 
and GOODS-South X-ray sources [heavy solid line), compared to the summed distribution of 
all the sources in model A [dot dashed line) and model B [light solid line). The distribution 
of unobscured (Type 1) [dashed line) and obscured (Type 2) [dotted line) AGN calculated 
using model B is also shown. The agreement is very good (K-S test gives 71% for model B 
and 53% for model A), showing that in a hard X-ray luminosity function both obscured and 
unobscured AGN are well represented. In the unified model (model B), the bright end of 
the distribution includes both Type 1 and Type 2 AGN, while the faint end is dominated by 
Type 2 AGN, as expected from the obscuration of the optical emission of the central engine. 
In all cases, most of the optical emission for Type 2 AGN comes from the host galaxy. 
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Fig. 8.— Redshift distributions for AGN in the GOODS-South (left panel) and North (right 

panel) fields. The observed redshift distribution {heavy solid line) includes both spectro- 
scopic {hatched area) and photometric redshifts and is 100% complete in the GOODS-S field 
and 75% complete in the GOODS-N region. The expected distribution {dashed line) for 
the complete AGN population detected in X-rays, calculated from our models (it is indis- 
tinguishable for models A and B), is similar but has more AGN at high redshift, especially 
compared to the distribution of spectroscopic redshifts. The distribution for sources with 
R < 2A mag {long dashed line) is very similar to the observed distribution of sources with 
spectroscopic redshifts. The discrepancy is larger in the GOODS-N field because of incom- 
pleteness in the redshifts. Spectroscopic redshifts are limited to AGN with R < 24 mag and 
thus exclude the high-redshift obscured AGN in the model, a larger fraction of which are 
included in the photometric-redshift sample. The data are therefore consistent with a signif- 
icant high-redshift population of obscured AGN which are missed in spectroscopic samples 
due to a selection effect. 
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Fig. 9. — Neutral hydrogen column density versus i-band magnitude for X-ray sources with red- 
shifts in the GOODS-South and -North fields. The Njj values were derived assuming an intrinsic 
power-law spectrum with photon index T = 1.9 and photoelectric absorption cross sections given by 
Morrison & McCammon (1983), shifted to the rest frame using spectroscopic redshifts from Szokoly 
et al. (2004) and Barger et al. (2003). The horizontal line at Nh = 10^^ cm~^ separates Type 1 
(unobscured) and Type 2 (obscured) AGN. The vertical line at i = 24.5 marks the spectroscopic 
limit for most cases. Given the spread in the intrinsic luminosity of the AGN and the constant 
contribution of an unreddened host galaxy, the correlation between Nh and i magnitude is weak. 
However, a general trend can be seen, with bright optical sources tending to be unobscured while 
those with higher values of A^^^ are fainter in the optical. Also, no Type 1 AGN is observed at 
magnitudes dimmer than i ~ 24, consistent with the predictions of our model (see Fig. 7). 
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Fig. 10. — X-ray-to-optical flux ratio distribution for sources in the GOODS-North and - 
South fields (solid hne) compared to the predicted distribution of model B (dashed line) . The 
X-ray flux is calculated in the (0.5-8 keV) band in the observed frame and is not corrected 
for obscuration. Optical flux is calculated from the z-band magnitude. Both distributions 
are very similar, except at the extremes. At low -Fx/-^opt non-AGN X-ray sources, mainly 
starburst galaxies, contribute to the observed distribution, while at the high Fx/F^^i, the 
most obscured AGN in the model are missing from the X-ray surveys. 
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Fig. IF — Predicted AGN number counts at 24/y,m for tlie total area and depth expected in the 
GOODS Spitzer Legacy fields (0.08 deg^). [long-short dashed line:) total AGN counts; [dashed 
line:) unobscured AGN; (dotted line:) obscured AGN; {solid line:) X-ray detected sources. For 
total and obscured counts, the lower line correspond to the infrared dust emission model with 
parameters Ri/Ro = 30 and q = I, while the upper line was obtained using Ri/Rg = 30 and q = 2. 
A third model with Ri/Ri = 100, q = 2 was also used and the results are between the other two 
models. For unobscured sources, our results are independent of the model used, which reflects 
the small dependence of the IR emission on the details of the torus composition and geometry for 
sources in which the line of sight does not intercept it. A significant number of the total AGN in 
the field, roughly 50% at the GOODS 24-/xm flux limit, are not detected by Chandra in X-rays. 
Most of these sources are obscured AGN with Nh > 10^^ cm~^. 
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Fig. 12. — Same as Figure 11, but for the S—fim IRAC band. In this case, for the total 
X-ray counts just one infrared model {Ri/Ro = 30, g = 1) was used since the difference with 
other parameters was small. Again, roughly 50% of the predicted Spitzer sources at the flux 
limit are not detected in the current Chandra deep fields. 
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Fig. 13. — Same as Figure 11, but for the 3.6-/im IRAC band. About 50% of the predicted 
Spitzer sources are missed in X-rays at the GOODS flux hmit. 



